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ABSTRACT 

We report Chandra observations of the complete sample of millisecond pulsars (MSPs) with precise radio 
positions in the globular clusters 47 Tuc (NGC 104) and NGC 6397. The x-ray luminosities and colors are derived 
or constrained and compared to x-ray MSPs previously detected in the field as well as one previously detected 
in a globular cluster (M28). The 47 Tuc MSPs are predominantly soft sources suggestive of thermal emission 
from small (r^ < 0.6km) polar caps on the neutron star rather than magnetospheric emission and are a relatively 
homogeneous sample, with most x-ray luminosities in a surprisingly narrow range (L, -1^ X 10^° ergs"'). We 
use previously derived intrinsic P values and find a new relation between L^ and spindown luminosity, E : Lj cx 
E ^, with (3 -0.5±0.2 vs. ~L0 for both pulsars and MSPs in the field. Adding the single MSP in NGC 6397 
constrains /3 = 0.5 ±0.15. This -E relation and also the Lj /E vs. spindown age are each similar to that 
found by Harding & Muslimov (2002) for thermal emission from polar cap heating. However, the cluster MSPs 
are relatively longer-lived (in thermal x-rays) than either the models or field MSPs, which may have additional 
magnetospheric (non-thermal) components. We suggest the cluster MSPs may have altered surface magnetic 
field topology (e.g. multipole) or their neutron stars are more massive from repeated accretion episodes due to 
encounters and repeated exchange interactions. MSP binary companions on or just off the main sequence (e.g. 
NGC 6397) are likely to have been re-exchanged and might show anomalous P and E values due to relaxation of 
misaligned core-crust spins. The radial distribution of ^40 soft Chandra sources in 47 Tuc is consistent with a 
~ L4 Mq component in a multi-mass King model and with the identified MSP sample. The implied total MSP 
population in 47 Tuc with Lj ^10^'^ erg s"' is ~35-90, and can constrain the relative beaming in radio vs. soft 
x-rays. NGC 6397 is relatively deficient in MSPs; its single detected example may have been re-exchanged out of 
the cluster core. 

Subject headings: globular clusters: general — globular clusters: individual (47 Tucanae) - 
x-ray: stars — binaries: general — pulsars: general 



stars: neutron - 



L INTRODUCTION 

Millisecond pulsars (MSPs) are end-products of binary evo- 
lution. In the Galactic disk, they are the fossil markers of rare 
primordial binaries in which neutron stars (NSs) are retained 
after birth in systems with low mass companions, from which 
accretion and a low mass x-ray binary (LMXB) phase trans- 
fers angular momentum to spin them up to msec periods. The 
"standard model" of MSPs as the remnants of LMXBs (re- 
viewed by Bhattacharya & van den Heuvel 1991) is supported 
by the discovery of the msec x-ray pulsar XTE J 1808-369 (Wi- 
jnands & van der Klis 1998) as well as the accumulating evi- 
dence for msec pulsation periods in a number of LMXBs from 
apparently stable frequencies in their quasi-periodic oscillation 
power spectra. However formation and evolution questions re- 
main since no LMXB in quiescence (qLMXB) has yet been 
found as a radio MSP, and the reduction in angular momentum 
transfer efficiency imposed by likely advection dominated ac- 
cretion flow (ADAF)-type accretion in quiescence means that 
the MSPs predominantly arise from a population of bright 
LMXB ti-ansients (Yi & Grindlay 1998) which may not be ob- 
served in numbers sufficient to support the disk MSP popula- 
tion. The recent radio detection of 20 MSPs in the globular 
cluster 47 Tuc (NGC 104) by Camilo et al. (2000), and the sub- 
sequent x-ray detection with Chandra (Grindlay et al. 2001a; 



hereafter GHEOla) of the 15 that were precisely located with 
radio timing solutions (Freire et al. 2001a; hereafter FCLOl), 
have opened a new window for the study of the nature and for- 
mation/evolution of MSPs. 

In globular clusters, in which MSPs are expected to show 
the same factor of 100-200 excess in their numbers (per unit 
stellar mass) as do the LMXBs compared with the disk pop- 
ulation (Clark 1975), the MSPs trace more exotic binary evo- 
lution (and possibly stellar evolution) paths than in the field. 
Here the MSP progenitor binaries can be formed readily from 
main sequence binaries by the exchange of a NS into the binary 
(e.g. Rasio, Pfahl, & Rappaport 2000) or they may also (in 
very high density cluster cores) form from direct 2-body tidal 
capture (Mardling 1995) or 3-body encounters. Clusters also 
allow the exchange of the original mass donor with a replace- 
ment main sequence (or near turnoff, given mass segregation) 
star, for which the MSP discovered in NGC 6397 (D'Amico 
et al. 2001a; D'Amico et al. 2001b; hereafter DPMOl; Ferraro 
et al. 2001) and detected with Chandra (Grindlay et al. 2001b; 
hereafter GHEOlb) might be the first example. 

The 15 MSPs with precise locations (FCLOl) include 47 Tuc- 
G and -I separated by just 0" 12. The sample is now extended to 
16 with the timing position for 47 Tuc-S (Freire 2001), which 
places it just 0''7 from 47 Tuc-F. We have identified (GHEOla) 
probable x-ray emission from the original 15 MSPs (-G and -I 
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cannot be resolved, and -C is a marginal identification at best, 
as described below), and we report here the possible detection 
of 47 Tuc-S. In this paper, we provide the details for the Chan- 
dra detections of this complete (for radio detection) sample of 
MSPs and derive the x-ray properties of this first significant 
sample of MSPs detected in a single globular. We derive the 
x-ray luminosities, spindown energy loss E, and the character- 
istic age {P/2P) for these MSPs, using the 47 Tuc distance, 
4.5±0.3 kpc, derived from the recent compilation of Zoccali et 
al. (2001). This is the first derivation of these properties for a 
significant sample of MSPs at a common known distance. We 
also compare the 47 Tuc MSP results to those of the single MSP 
detected in NGC 6397 (GHEOlb) and previous (ROSAT) re- 
sults for MSPs in the field and one other globular (M28), as 
summarized in the review by Becker & Triimper (1999; here- 
after BT99). 

2. CHANDRA OBSERVATIONS AND ANALYSIS OF MSPS IN 

47 TUC 

A 74ksec Chandra observation (72ksec total data obtained) 
was conducted on 2000 March 16-17, in which at least 108 
sources were detected in the 2' x 2. '5 around the cluster cen- 
ter (GHEOla). With the ACIS-I detector used, the 0."492 pixel 
size and 0."7 (FWHM) telescope-detector point spread function 
(psf) for sources within ~4' of on-axis combine to give source 
positions with la centroid uncertainties of <,0"3 for sources 
detected with >,90% confidence (Jerius et al. 2000). Even a 
source identified with just 1 count (cf. discussion below of 
MSP 47 Tuc-C) is thus located with ~ 0."5 (Ict) centroiding 
error. The uncertainties in the absolute positions of Chandra 
sources are dominated by aspect (guide star) offsets (i.e. bore- 
sight offsets between the x-ray and optical axes) and are cur- 
rently estimated to be ~ Q'.'l (Aldcroft et al. 2000). 
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Fig. 1 . — Images (0.5^.5 keV) of the MSPs in 47 Tuc not detected (or 
confused) by WAV detect (Q, C, T, G & I, F & S, and L), showing probable 
emission from the MSP in each case, and the weakest detection (M). Pixel ran- 
domization has been removed from the pipeline processing, and 1 " circles are 
placed on the radio MSP locations. These circles are used to extract the counts 
for each MSP, as they contain ~ 95% of the encircled energy for ^0.3 keV 
sources, which describe all of the MSP spectra except for -J (cf. Fig. 2). 

Given these relative and absolute errors for source positions, 
the initial run of the CIAO"' tool WAVDETECT (Dobrzycki et 
al. 1999) in our general study (GHEOla) of the sources near 
the center of 47 Tuc detected 10 sources containing 11 of the 

All CXC software is available at http://asc.harvard.edu/ciao. 



MSPs with ~3-10mas positions as given by FCLOl. A con- 
sistent offset of (RA, Dec) = (0"38±0"14,-0."82±0."11) was 
found for 6 relatively isolated MSPs (47 Tuc-D, -E, -H, -J, -N, 
-U) detected with wavdetect. Applying this offset to correct 
the Chandra positions to the MSP reference frame, we derived 
a small linear correction of the Chandra coordinates in both RA 
and Dec to minimize the rms deviations (with RA, Dec values 
0"08,0f'll for the positions given in Table 1) with the radio 
positions. Details of this astrometry are given by Edmonds et 
al. (2001), who extended the astrometric solution to HST data 
to identify the hkely optical counterpart for 47 Tuc-U. For the 9 
remaining sources of FCLOl, including the unresolved pair -G 
and -I, we used the precise astrometric solution to examine the 
Chandra images (in each band) for evidence of these MSPs (cf. 
Fig. 1). For 3 of these (47 Tuc-L, -Q, -T), WAVDETECT failed to 
detect the source due to crowding with a brighter neighbor, al- 
though inspection showed clumps of photons (similar in color 
and flux to the detected MSPs) centered at the radio position 
and consistent with the expected psf. For 47 Tuc-C only one 
count is detected at the predicted MSP position, but the rela- 
tively large exposure correction (factor of 2.5 [cf. Table 1], due 
to its being near a detector gap) for this source is consistent with 
it having x-ray flux within a factor of ^2 of the faintest other 
MSP, 47 Tuc-M. 

Total counts were estimated for each MSP by extracting 
counts in each band in a weighted (by the psf) 1 " radius cell 
centered on the MSP position and allowing for psf overlap 
contributions from neighboring sources using their WAVDE- 
TECT counts in each band and offset positions. The recently 
derived position for 47 Tuc-S (Freire 2001) shows that the orig- 
inal WAVDETECT flux (and thus derived x-ray colors and L^i) 
for 47 Tuc-F (Table 1 of GHEOla) is contaminated by 47 Tuc- 
S, separated by just 0"7. wavdetect identified an elongated 
source between the positions of 47 Tuc-F and 47 Tuc-S, which 
can be best explained as the combined emission of the two 
sources. We estimate the counts for -F and -S by noting that the 
x-ray position is twice as close to -F as -S, and thus assign -F 
two-thirds of the counts. This is supported by visual inspection 
of the images in each band. Detected position offsets in the 0.5- 
4.5 keV primary detection band for the 10 wavdetect sources 
(12 MSPs) and counts in 3 bands (see below) for all 16 MSPs, 
are given in Table 1. 

2. 1 . X-ray Colors and Emission Models 

Our initial analysis of the MSP sample indicated the objects 
are relatively soft sources compared with the probable cata- 
clysmic variables (CVs; GHEOla). In Table 1 we give the de- 
tected counts in 3 bands: softcts (0.2-1 keV), mediumcts (1- 
2keV) and hardcts (2-8 keV) for each of the 14 resolved MSPs, 
with counts for 47 Tuc-G and -I (unresolved) estimated. We use 
3 bands rather than the hardness ratio Xcolor = 2.51og(cts[0.5- 
1.5keV]/cts[1.5-6.0keV]) used in GHEOla since it allows an 
approximate spectral analysis for the MSPs. With (limited) 
counts in 3 bands, we form the hardness ratios HRl = medi- 
umcts/softcts and HR2 = hardcts/mediumcts and plot the MSPs, 
with counting statistics errors, in the color-color diagram shown 
in Figure 2. The MSP colors are clustered in a relatively narrow 
range of HRl and HR2; 47 Tuc- J is clearly (much) harder, as 
was evident in the Xcolor distributions plotted in GHEOla. 

For comparison with the data, we construct (using the 
PIMMS tool for ACIS-I) values of HRl and HR2 for 3 sim- 
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pie models derived for the 47 Tuc absorption column (Nh = 
2.4 X 10^° cm-2; cf. GHEOla): thermal bremsstrahlung (TB), 
blackbody (BB) and power law (PL), with index values (kT or 
photon index) given in the caption of Figure 2. It can be seen 
that the observed range of HRl vs. HR2 is roughly consistent 
with TB spectra with kT ~ 1 keV, BB spectra with kT ~0.2- 
0.3 keV or PL spectra with photon index -^3. The track for BB 
models is most consistent with the data for all but 47 Tuc-J. 
The weighted mean colors are plotted for all the MSPs and for 
all but 47 Tuc-J. This is consistent with a single kT ^0.22 keV 
for all the MSPs (except -J). Thus we give in Table 1 approx- 
imate x-ray luminosity values for a BB spectrum with kT = 
0.22 keV. Apart from 47 Tuc-C, all other MSPs are detected 
with in a surprisingly narrow range: ^1-4 xlO^" erg s"', 
where the Lj values are again (as in GHEOla) quoted for the 
0.5-2.5 keV band to facihtate comparison with ROSAT results 
for MSPs (BT99) and to minimize uncertainties due to extrap- 
olation of uncertain spectra over a broad band in which counts 
are not actually detected. 
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Fig. 2. — X-ray color-color diagram for all located MSPs (except 47 Tuc-C 
which is too faint, and -G which is combined with -I) in 47 Tuc (solid dia- 
monds) vs. tracks for models with parameters changing from lower left to up- 
per right: blackbody (BB; open squares), with kT = 0.2, 0.25, 0.3 keV; power 
law (PL; open triangles), with photon index a = 3, 2, 1.5, 1; and thermal 
bremsstrahlung (TB; open circles), with kT = 1, 2, 3, 6keV. Weighted mean 
colors for the MSPs in 47 Tuc with and without 47 Tuc-J (labeled separately) 
are shown as the solid square and triangle, respectively. 

Several constraints on MSP x-ray emission may be inferred 
from Figure 2. Given the average MSP luminosity (from Table 
1) of Lj: ~2 X 10^" erg s~^ then for a BB with the suggested 
kT ~0.22keV or T ^2.6 x 10''K, the bolometric luminosity is 
i^x-boi ~3 xlO^" erg s"'and the total emitting radius is only 
Tx ~0.2km. NS H-atmosphere models, as needed for thermal 
spectra of the quiescent LMXBs in 47 Tuc (cf. Heinke et al. 
2002 and references therein), would decrease the temperature 
by a factor of ^2 and thus increase the emission radius to ^0.8 
km. For a dipole field, both thermal (isotropic) polar caps are 
likely visible and the effective radius of each would be ~0.14 
km, or '~0.56km for H-atmosphere models. Such a small emis- 
sion area strongly suggests emission from a radiating area much 
smaller than the entire polar cap, which has expected radius of 
Tpc = (fiR*/c)'/^R* for a MSP with angular rotation frequency 
O and NS radius R* (or tpc = 2.6km for a 3 msec period MSP 
and R* = 10km). Partially filled polar cap thermal emission (i.e. 
Ta < Tpc) was also found by Zavlin et al. (2002) for the bright 
MSP J0437-4715 (for which a two-component model could be 
fit to the much better statistics data than possible for the 47 
Tuc MSPs) and is predicted by recent polar cap heating models 



of Harding & Mushmov (2002; hereafter HM02), although we 
find important differences (cf. §5). The small emission area 
implies x-ray pulsations (which cannot be measured given the 
3.2 sec ACIS-I integration time) with a sinusoidal pulse shape 
appropriate to the fractional visibility of the isotropically radi- 
ating thermal polar cap(s). In contrast, the narrower pulse duty 
cycles of ~10% for some field MSPs (and one in the globu- 
lar cluster M28; cf. BT99) are probably due to non-thermal 
beamed emission. 

If the x-ray emission were a (pure) TB spectrum with kT 
~ 1 keV, the required emission measure EM ~4 xlO^^cm-^ 
would suggest plasma densities n ~3 x lO^^R"^/^ cm~^ for an 
emission region (assumed spherical) of radius R. Such a region 
must not dominate the observed dispersion measure of the MSP, 
which is typically DM = 24.4 cm"^pc for the 47 Tuc MSPs 
with a variation SDM ~0.1cm"^pc that is likely due to the 
ionized ISM in 47 Tuc (Freire et al. 2001b). Thus, the MSP 
local contribution to DM is probably n • R ^0.03 cm"^pc ^1 
X 10'^ cm"^. Combined with the constraint above for EM, for a 
constant density emission region this impUes an unreaUstically 
large lower bound for the emission region: R ^3pc. (For a TB 

source with radius R ~10^" cm, or comparable to the size of 
the typical MSP binary separation, n -^3 x lO'^ cm"^ is re- 
quired which would strongly violate the DM constraint.) Thus 
a pure TB spectrum is ruled out. Similarly, with the exception 
of 47 Tuc-J, PL models are unlikely since their photon index 
must be implausibly large (>,3) and are still inconsistent with 
the colors for most of the sample. The limited statistics can not, 
of course, rule out combined BB -i- PL models, but certainly it 
appears that the MSP sample in 47 Tuc is dominated by thermal 
emission from polar caps of the NSs. 

3. X-RAY vs. RADIO PROPERTIES OF MSPS IN 47 TUC 

The 47 Tuc MSPs are the first significant sample of MSPs 
at known distance with measured x-ray and radio properties. 
Several long-standing pulsar questions can be addressed. 

3.1. X-ray vs. Pulsar Spin Properties 

Foremost is the correlation of x-ray luminosity and pulsar 
spindown luminosity E, which is found for field MSPs (with 
much more uncertain distances) to scale as L^ ^lO""' E (BT99). 
Here the globular cluster is both a hindrance and an interesting 
laboratory: the MSPs are accelerated in the cluster potential, 
thus biasing the observed values of P and hence E (cx P/P^). 
FCLOl have modeled the variations of P among the 15 MSPs 
with precise locations (not including 47 Tuc-S) to derive con- 
straints on both the cluster potential and the intrinsic val- 
ues for the MSPs. Tighter constraints on P, can be obtained 
by assuming the cluster gas (detected from variations in DM) 
is uniform so that the observed DM value and projected radius 
give the 3D position in the cluster for each MSP (see Freire 
et al. 2001b). Assuming a King model for the cluster, and us- 
ing the central velocity dispersion (7^0= 11.6kms~' (Meylan & 
Mayor 1986), then gives the cluster acceleration term (cf. Fig- 
ure 2 of Freire et al. 2001b) which is subtracted from the ob- 
served P/P to infer the intrinsic P,. Using a standard NS mo- 
ment of inertia / = 10^^ g cm^, we then derive E = Air^IPi/P^ 
for each MSP. Results are given in Table 1 and plotted vs. L^ in 
Figure 3. Uncertainties in the derived E values are typically 
0.2-0.5 in the log (though several are larger) but are not shown 
for clarity; uncertainties in log(L;c) are typically 0.2 (see Table 
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1). For comparison with all MSPs previously detected in x- 
rays, we also plot values for 10 field MSPs and one each in the 
globular clusters M28 and NGC 6397. The data for the field 
MSPs and M28 (from which E values are derived) are summa- 
rized in the new compilation of data and distances in Table 2 
whereas the NGC 6397 x-ray results are from GHEOlb and are 
discussed in detail in §4 below. We have repeated this analysis 
for variations in both cluster distance (S.Okpc; +\.5a) and cen- 
tral velocity dispersion (13.0km s"'; +\a) to sample the range 
of cluster acceleration models, and thus intrinsic P, values. Our 
goal is to derive the approximate index (3 in the relation hx cx 
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Fig. 3. — Lv VS. E for MSPs in 47 Tuc (diamonds) and field (triangles). 
MSP in NGC 6397 (square) is labeled as is -J in 47 Tuc. Values for MSPs in 
field and M28 are from Table 2, with = 4ncPF, and E = 4n^IPi/P^. 

Whereas the field MSPs (and pulsars generally) are roughly 
consistent with the linear relation ^10"^ E (BT99), the 
MSPs in 47 Tuc appear to have a much weaker dependence: 
logU = (0.47 ± 0.10)log£' + 14.2 ±3.4 for the nominal cluster 
model with cr,,o= 11.6kms"', where the errors (±1<t) are due 
to just the scatter in the points. Including the errors on each 
point, but with the approximation that unequal errors (on E) are 
simply averaged (which biases the slope to steeper values, since 
the unequal errors are much larger for smaller values of E), in- 
creases the slope to 0.59 ±0.28 and offset to 10.0 ± 10.6. The 
best estimate of the slope is thus /3 ~0.5 ±0.2. Allowing for the 
error in the cluster model by using the -i- 1 cr value for the cen- 
tral velocity dispersion, (t,,o= 13.0kms"', and again including 
the errors on each point gives slope 0.48 ±0.21 and intercept 
13.8 ± 7.5, and so is consistent. Similar variations in p and the 
intercept are found for the increased (by about 1 .5cr) cluster dis- 
tance, 5kpc. Thus, the MSPs in 47 Tuc have L, (0.5-2.5 keV) 
values within a factor of ^4 (apart from 47 Tuc-C) despite a 
range of ~25 in E. For comparison, the field (and M28) MSPs 
in Figure 3 have (3 = 1.13±0.15 and offset -7.48 ±5.03 and 
so have /3 consistent with both the ROSAT band (0. 1-2.4 keV) 
results of BT99 and the 2-10 keV results derived by Possenti et 
al. (2002), who find f3 = 1.37±0.10. 

It is tempting to speculate these different f3 values may be 
due to the different formation histories and possibly different 
physical parameters of MSPs in the field and 47 Tuc. In Fig- 
ure 4a we plot vs. spindown ages, r = P/2P,, for field and 
47 Tuc MSPs. Error bars on the age parameter are not shown, 
for clarity, but are typically ^^JJ [ in the log and primarily due 
to the cluster acceleration model. Despite the uncertainties, the 
correlation is striking: the field MSPs show a declining with 



"age," and the 47 Tuc MSPs appear to fall on this trend but 
with a flatter L, vs. age slope. (Although both age and E are 
derived from combinations of P and P,, they are algebraically 
different so that the same flatter correlations with are not 
inevitable.) In §5 we compare the IE -r relation (Fig. 4b) 
with the models of HM02, finding both qualitative similarities 
(noted by HM02) but also important differences. 
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Fig. 4.— Top (a): L,- vs. spindown age (t = P/2P,) for MSPs in 47 Tuc 
(diamonds) and field (triangles). MSP in NGC 6397 (square) is labeled as is 
-J in 47 Tuc. Values for MSPs in field and M28 are from Table 2. Bottom (b): 
Lv /E vs. spindown age, and approximate correction for pulse period P (from 
HM02; see text for empirical dependences on P and t) for 47 Tuc vs. field 
(and NGC 6397) MSPs. 

We note that although spindown ages correspond only ap- 
proximately to actual ages, the derived values are consistent 
with the formation of most 47 Tuc MSPs in a "burst" early in 
the cluster history. However this by itself would not provide 
an explanation for the different correlations observed for Li in 
47 Tuc and field MSPs. As a measure of the systematic uncer- 
tainties in age, the 0.6 Gyr value derived from the cooling age 
for the He- WD companion of 47 Tuc-U (Edmonds et al. 2001) 
may be compared to the ^0.9 Gyr or 2 Gyr spindown age of the 
MSP (two solutions for P, are possible for 47 Tuc-U, as is also 
the case for -G and -M from the total sample of 15 for which the 
cluster acceleration model yields values for P, and thus E and 
r). 

Another possible physical difference between the 47 Tuc and 
field MSPs might be magnetic field strength at the NS sur- 
face and thus also the light cylinder (at which the corotation 
speed equals c), since this likely affects the relative importance 
of non-thermal (magnetospheric) emission. For an assumed 
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dipole field, this is given by = 9.35 x 10^P/''^Pmsec G, with P, 
in units of 10~^° s s"' . In Figure 5a we plot for the same pulsars 
vs. Bic- As before, considered as a homogeneous group, the 
field pulsars (and M28) lie on a significantly steeper logarith- 
mic slope (1.70±0.25) than the 47 Tuc MSPs (0.66±0.13). In 
contrast the field at the NS surface, B^nd = 3.2 x lO'*^ (PP)'/^ G, 
is nearly independent of . The correlation of Lj with Bjurf 
(Fig. 5b) is less defined, with twice the scatter for both field 
and 47 Tuc MSPs, and with logarithmic slopes differing even 
more: 0.05±0.27 for the 47 Tuc MSPs vs. 2.80±0.99 for the 
field MSPs. Whereas at the light cylinder the dipole component 
of the field, Bic, would dominate, the surface field would be 
sensitive to multipole components which could give the larger 
scatter in Figure 5b. We return in §5 to additional arguments 
for a modified surface field topology. 



1821-24/M28 



Log(Lx) = 1.7aLog(BJc) + 23.0 (Field & W2i) 



Log(Lx) = 0.66Log(B Jc) + 27.2 {47Tuc) 



5.00 5.50 
Log(B_lc, gauss) 



Log(Lx) - 2.SaLag 



Li}9(U) = 0.05Log(B_surf) t 29.8 (47Tuc) 



7.50 8.00 8.50 9.00 9 50 10.00 10 50 

Log(B_surf, gauss) 

Fig. 5. — Top (a): Lv vs. B[c, inferred magnetic field strength at tlie liglit 
cylinder (ric = cP/2n) for MSPs in 47 Tuc (diamonds) and field (triangles). 
MSP in NGC 6397 (square) is labeled as is -J in 47 Tuc. Values for MSPs 
in field and M28 are from Table 2, with Bi^ derived assuming standard 1 / P 
dependence for a dipole field (see text). Bottom (b): Lv vs. surface magnetic 
field (see text). 

We note that 3 out of the 4 MSPs in the field with Bi^ > 
10^'^ G display x-ray emission that is seemingly magneto- 
spheric, with the nature of emission from the 4th (the eclipsing 
MSP B1957H-20) indeterminate (BT99; Takahashi et al. 2001). 
Conversely, field pulsars with Bi^ < 10^ G have x-ray emission 
that is typically either thermal or of indeterminate character 
(BT99). Considering the small numbers of such pulsars stud- 
ied, and that most of them have highly uncertain distances (see 
Table 2), it is possible that field pulsars with Bic < 10^ G may 
show an Lj-Bk trend that is fairly flat and roughly consistent 
with the better determined relation for the 47 Tuc MSPs. How- 
ever for this interpretation to hold, a few field MSPs with well 



determined distances (e.g. J0437-4715 and J1744-1134; cf. 
Table 2) must be accounted for and the even larger deviations 
of J0751H-1807 and J 1024-07 19 from the 47 Tuc coiTelation 
line would require a factor ^3 adjustments in these MSP dis- 
tances. The fact that 5 of the 7 field MSPs with log(£') <34 
have Lf values above the trend-line for the 47 Tuc MSPs may 
suggest these field MSPs still have non-thermal emission com- 
ponents. We return (in §§5 and 6) to further discussion of the 
possibility that the magnetic field configuration and evolution 
for old cluster NSs may differ from those in the field. 

3.2. Radial Distributions 

FCLOl pointed out the puzzling fact that the radial distribu- 
tion of the (then) 15 MSPs with timing positions seemed to be 
truncated within '^3.5rc, where = 24" is the core radius of 
the visible stars in the cluster as determined most recently (and 
completely) by Howell, Guhathakurta, & Gilliland (2000). The 
addition of 47 Tuc-S, as the 16th MSP with a precise location 
(Freire 2001), does not change this. It is therefore of partic- 
ular interest to compare the radial distribution of the 16 MSPs 
with the Chandra source distribution over a region large enough 
to measure the radial profile of the sources most likely to be 
MSPs. Our 2' X 2'.5 centi-al field study (GHEOla) suggested 
most of the MSPs are soft and a possible total MSP popula- 
tion of ^^50 in this region alone. Here we extend the study of 
the Chandra sources out to a radius of 4' (^10 r^.), to search for 
possible cutoffs or truncations in any of the Chandra source dis- 
tributions as might be expected if the MSP distribution were cut 
off. We compare the MSP radial profiles with those of Chan- 
dra sources in several apparent spectral (hardness) classes to 
further test the possible identification of faint soft sources ("red 
background sources") with the large MSP population suspected 
(Camilo et al. 2000; GHEOla). This analysis will include, then, 
the vast majority of sources associated with 47 Tuc. Complete 
details (spectra, timing) of this source population will be in- 
cluded in the still-larger study (fuH ACIS 16' field) of 47 Tuc 
to be reported by Heinke et al. (in preparation). 

WAVDETECT analysis of the the 4' (radius) field yields a to- 
tal of 187 sources above a limit of medcts(0.5^.5keV) = 3 
(vs. 103 in GHEOla within the central 2' x 2.'5 box), or 126 
sources above a 5ct limit (vs. 92 in the 2' x 2'.5 box), which we 
impose to get a > 3(t significance. For softness ratio (Xcolor 
= 2.51og[Soft(0.5-1.5keV)/Hard(1.5-6.0keV)]; cf. GHEOla) 
determination, we impose the WAVDETECT count limit of 3cts 
as an upper limit in either Soft or Hard bands if there is no 
WAVDETECT detection in that band. This alters the Xcolor val- 
ues somewhat (for weak sources) from those given in GHEOla. 
In Figure 6a we provide the radial profiles of the sources de- 
tected with x;,5cts in at least one band, which are divided into 
Soft (Xcolor ^1.0; 44 sources) and Medium (-1.0 ^Xcolor 
^1.0; 41 sources) groups. For the Medium group we also re- 
quired that the source be detected by WAVDETECT (^3 cts) in 
both the soft (0.5-1.5 keV) and hai'd (1.5-6.0keV) bands. Not 
plotted are the Hard (Xcolor ^-1.0; 9 sources) sources and re- 
maining sources (32) not included in the above 3 groups, since 
the Hard sources likely have a significant background AGN 
contribution (though at least 2 of these 9 sources are cluster 
CVs; GHEOla), and the remaining sources are indeterminate. 
The Medium source group, and thus its radial distribution, is 
also likely to be contaminated in this larger field by background 
AGNs. Scaling from the deep survey counts and AGN colors 
measured by Giacconi et al. (2000), we expect ~8 Medium 
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Fig. 6. — Top (a): radial source distributions for 47 Tuc MSPs vs. Chandra Soft (not including MSPs) and Medium (raw and with statistical subtraction of expected 
AGNs) sources detected with >,5cts together with best-fit King model predictions for each group. Bottom (b): envelope (5% and 95%) of radial distributions for 
1000 samplings of 16 soft Chandra sources (hatched region) vs. the MSP sample (soUd histogram). 



and ^3 Hard sources to be background AGN. Therefore, the 
Medium source distribution has also (Fig. 6a) been corrected by 
the following Monte-Carlo procedure. A total of 1000 bootstrap 
resamplings of the distribution were generated. For each of 
these samples, an expected number of background sources was 
selected from a Poisson distribution with mean 8. This number 
of radial source positions was then selected from a uniform spa- 
tial distribution (i.e. with cumulative distribution oc r^), and the 
closest actual sources to these positions were removed from the 
sample. The average of the 1000 background-corrected sam- 
ple distributions was then adopted as the best estimate for the 
corrected cumulative distribution. 

In Table 3 we give the probabilities from a KS test that each 
source group has a radial distribution consistent with that of 
any other group (including the 16-memberMSP sample). For 
the Medium group, we show results for both the raw and cor- 
rected distributions. All of the cross-sample comparisons indi- 
cate consistency (i.e. probability of the samples being consis- 
tent exceeds 5%) with the exception of the Medium-raw - Soft 
case. We ascribe this inconsistency to the contamination of the 
Medium-raw sample by background sources and note that the 
Medium-corrected sample is consistent with the Soft sample. 

In Figure 6b we show that 1000 repeated samplings of 16 
sources from the Soft population gives a range of radial distri- 
butions that includes the observed MSP distribution. We note 
that both the CVs and BY Draconis (BY Dra) main sequence 
binaries have radial distributions that are also consistent with 
that of the Soft sources. However the radial incompleteness of 
our present HST sample and thus CV and BY Dra identifica- 
tions makes this result preliminary. 

For each of the observed or corrected (Medium) radial distri- 
butions in Figure 6a, we have fitted what we shall refer to as a 
"generalized King model" profile (cf. Lugger, Cohn, & Grind- 
lay 1995) to the projected surface density, 



The corresponding space density profile is given by 
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Here a is the large-r power-law index and the core radius is 
related to the radial-scale parameter ro by 
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Following the analyses of Jernigan & Clark (1979), and 
Lightman, Hertz, & Grindlay (1980), we assume that the pro- 
jected distribution of the dominant visible stellar population, 
with mass (~ 0.7 Mq in the cluster core due to mass seg- 
regation), is described by an "analytic King model," i.e. the 
form given by eq. (1) with a* = -2. In thermal equihbrium, the 
space density of an x-ray source population with stellar mass 
= qM^ is given by ^(r) oc n*(r)^. Thus, the power-law in- 
dex of the projected x-ray source density profile is predicted to 
be 



and the corresponding core radius is predicted to be 



0.68 q 



-1/2 



(4) 

(5a) 
(5b) 



The latter approximation holds in the limit q^\ and is accu- 
rate to 7% for q>2. 

In Table 4 we list the values derived for r^. and a for each 
source group. The Soft sample is of particular interest, since 
it is the largest one (A^ = 44) and its parameters are reason- 
ably well determined. As we have shown in Table 3, both the 
Medium-corrected and the MSP samples are consistent with the 
Soft sample. In principle, the mean x-ray source mass may 
be computed from eqs. (4) and (5). Eq. (4) gives q=1.5 for the 
best-fit value of ax = -3.6 ± 0.8, with a l-a uncertainty range 
of ^ = 1.3-1.8. Eq. (5a) gives q=\.7 for the best-fit value of 
r„ = 15. "2 ± 3"5 together with r„ = 24"0 (Howell et al. 2000), 
with a 1-cr uncertainty range of ^ = 1.3-2.5. Thus, both mea- 
sures of q put the best-fit value at around 1 .6 and allow values as 
large as <7 « 2. Given our estimate that the dominant mass in the 
core of 47 Tuc is ^ 0.7 Mq, this implies that Mj lies in the 
range ~ 1.1-1.4M0. For comparison, Grindlay et al. (1984) 
obtained a most probable value of ^ = 2.6 for an ensemble of 
eight luminous LMXBs in as many clusters. 
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There are a number of limitations to the "classic" analysis 
outlined above, which may well affect the estimate of the x- 
ray source mass. First, this analysis assumes that the poten- 
tial is dominated by stars of a single mass. In fact, there is a 
large range of stellar mass with substantial mass segregation. 
Secondly, the analysis assumes that the x-ray sources make a 
negUgible contribution to the gravitational potential. However, 
given mass segregation, even a modest NS or heavy white dwarf 
(WD) population can make a substantial contribution to the po- 
tential in the innermost part of the cluster. Finally, the analysis 
assumes that the cluster is in a state of thermal equilibrium, i.e. 
that the cluster is isothermal and that each species is in energy 
equipartition. This condition does not hold rigorously, even in a 
nearly isothermal multi-mass King model. Also, the relatively 
low Mj ~1.1-1.4M0 may suggest non-equilibrium (perhaps 
from binary scattering) if the NS masses are >, 1.4 Mq as pos- 
sibly suggested by the long-Uved thermal emission (cf. §5). 

In light of these considerations, it is useful to compare our re- 
sults for the Soft source distribution with the predictions from 
actual dynamical models. Meylan (1988) fitted multi-mass, 
anisotropic King-Mitchie models to surface brightness and stel- 
lar velocities for 47 Tuc, in order to constrain cluster parameters 
including the mass function. The models include a NS com- 
ponent, which comprises between 0.1-0.5% of the total clus- 
ter mass of 7 X 10^ Mq. Radial profiles for all components are 
given for one model with a NS mass fraction of 0.2% (equiva- 
lent to ~ 1000 NSs in the cluster). The power-law for the NS 
component has a slope of -4.0 in projection, which is consistent 
with the value of a = -3.6 ± 0.8 found here for the Soft source 
distribution. The predicted core radius for the NS component 
is 16", which is also in good agreement with our determination 
of Tc = 15"2±3."5 for the Soft sources. By comparing the ob- 
served Soft source distribution with a set of models with a range 
of NS fraction, it should be possible to constrain this parame- 
ter. The models may be generated either by following Meylan's 
(1988) King-Mitchie model approach, or by Fokker-Planck in- 
tegrations (e.g. Murphy et al. 1998). We are undertaking a study 
based on the latter method. 

3.3. Candidate MSP Sample in 47 Tuc and Estimated Total 

The KS probabilities (Table 3) and King model fits (Table 4) 
show that the radial distribution of the 44 Soft sources is in- 
distinguishable from that of the MSPs, and indeed 7 of the 9 
sources identified with MSPs with well-measured Xcolor val- 
ues are soft. These soft MSPs include the doubles (-G, -I) and 
(-F, -S), for which separate Xcolors could not be determined, 
so that in fact 9 MSPs are accounted for. The remaining MSPs 
are either faint (-M, -T; softcts <,5; cf. Table 1) or crowded 
(-L, -Q), but it can nevertheless be seen from the x-ray colors 
plotted in Figure 2 that all 47 Tuc MSPs, with the exception of 
47 Tuc-J, are probably Soft (although 47 Tuc-U, with Xcolor = 
0.9, is formally in the medium group, while colors for 47 Tuc-C 
could not be measured). 

Motivated by these results we have searched the list of Soft 
sources for extra MSP candidates. Of the 37 Soft sources 
not identified as MSPs, two are likely qLMXBs (X5 and X7; 
GHEOla and Heinke et al. 2002), two are CVs (V3 and AKO 
9; GHEOla) and three are sources with counts more than two 
times brighter than the brightest detected MSP. Only 17 of the 
remaining 30 sources fall in the HST field of view (FOV) of 
Albrow et al. (2001). Nine of these 17 sources are astromet- 
rically matched to optical variables (hereafter "active binary" 



candidates) identified with HST (Edmonds et al., in prepara- 
tion), and an additional 2 are faint blue stars and possible CVs, 
one of which is also a periodic variable. Seven of the optical 
variables have independently been identified as likely BY Dras, 
red stragglers, or contact/semi-detached binaries by Albrow et 
al. (2001), and 2 of them were reported by GHEOla as BY Dras 
based on the additional (more restrictive) constraint of signifi- 
cant x-ray variability. 

Despite the lack of detection of any of the 7 soft MSPs as 
optical variables, it is possible that several of the active binary 
candidates contain MSPs. For example the optical companion 
to the one MSP known in NGC 6397 was previously identified 
as a BY Dra candidate by Taylor et al. (2001; hereafter TGEOl). 
The corresponding Chandra source is relatively hard (see § 4) 
and would not have been included in our list of soft sources in 
47 Tuc, but we cannot be sure that all other MSPs with main 
sequence (or red straggler, as likely for the NGC 6397 MSP) 
companions will be similarly hard. In fact, after the initial sub- 
mission of this paper, the relatively hard Chandra source W29 
(GHEOlb) has been identified with the MSP 47 Tuc-W as hav- 
ing a main sequence companion (Edmonds et al. 2002; cf. § 4). 
Also, one or both of the soft sources with a faint blue opti- 
cal counterpart could be a He WD companion to an MSP (like 
47 Tuc-U; Edmonds et al. 2001), instead of a CV. The definitive 
MSP characteristic is, of course, a pulsed radio source. How- 
ever, a significantly larger sample of MSPs with radio timing 
positions, for comparison with the optical variables, will be 
difficult to obtain because so much observing time has been de- 
voted to finding the 16 already known. Therefore, to estimate 
the number of MSPs implied by our soft sample we consider 
two extreme cases: (a) none of the candidate optical counter- 
parts (9 active binary and 2 CV candidates) are MSPs, or (b) 
all of the candidate optical counterparts are MSPs. For case (a) 
we find that only 6 of the soft sources falhng in the HST FOV 
(and with counts < 2 x the brightest detected MSP) do not have 
a known ID. Correcting for the detected fraction of MSPs with 
timing positions that are soft sources (7/16) and the incomplete 
FOV of HST (17/44 sources are outside the HST field) we es- 
timate an extra population of 19 (= 6[1 + 17/44] x 16/7) MSPs 
with L;f > 10^^ erg s ', in addition to the 16 already known. 
For case (b) a total population of 69 (= 30 x 16/7) extra MSPs 
is imphed, for a total maximum MSP population of 85. 

These estimates are based on scaUng from the number of soft 
source detections with wavdetect to the total MSP popula- 
tion with timing positions. Using a different technique we esti- 
mate the number of soft sources missed by WAVDETECT by ex- 
amining the net counts not attributable to wavdetect sources 
within a radius of 24" (=1 for the background stars), where 
crowding will most limit faint source and crowded-field detec- 
tion, and the annulus 1-2 surrounding the cluster center. By 
removing counts from the original image in 1 " radius regions 
around each wavdetect source (from GHEOla) in both the 
Softcts and Mediumcts bands (cf. Table 1), and then subtract- 
ing a fraction of the remaining counts (in each band) to correct 
for the psf-spillover from these bright(er) sources due to the 
energy-dependent Chandra telescope psf, we find 180±56 Soft- 
cts, 20±58 Mediumcts and 15±36Hardcts in the central 1 ra- 
dius region, and 60±24, 75±28 and 45±25, respectively, in the 
surrounding 1-2 annulus. Thus, as suspected in GHEOla, 
there is an excess of Softcts in the cluster core, but the uncer- 
tainties in this extraction are large and no excess in the core of 
Mediumcts or Hardcts is required. For the detected MSPs with 
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Softcts as listed in Table 1 (with mean value 6cts), the excess 
total Softcts (240) could arise from an additional ^40 MSPs, 
all with Lj: ^,10^*^ erg s~^ Since some of these soft sources 
could be BY Dras (coronal sources are typically soft) or CVs, 
as already noted (GHEOla), we regard 40 as an upper limit and 
estimate ~20 (half) "unresolved" MSPs in the core. 

Scaling by our known sample of 16 MSPs, we derive 40 x 
(16/7) = 91 MSPs as an upper limit on the total MSPs in the 
cluster. This upper limit agrees with the upper limit derived 
from consideration of WAVDETECT soft sources, and we adopt 
90 as an approximate maximum number of MSPs in 47 Tuc 
with L;c ^.lO^" erg s"'. Both of these estimates assume that no 
beaming of the MSPs in x-rays is occurring. (Note that thermal 
polar cap emission, while pulsed, is not beamed.) If there are 
more MSPs than this in 47 Tuc, they must have relatively hard 
spectra (since the above totals are scaled from the soft source 
counts) which could make them like the MSP in NGC 6397 
(see below) for which radio emission is occulted much of the 
time. Alternatively, additional MSPs are beamed away from us 
in the x-ray, or of significantly lower x-ray luminosity on av- 
erage. The softer color of the excess counts than our known 
MSP sample suggests that either the excess counts are com- 
posed largely of other sources (such as faint BY Dra binaries), 
or that there is a significant population of faint (in x-ray and 
radio) MSPs which would have cooler (and thus softer) polar 
caps. 
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Fig. 7. — hx VS. Lra(iio for MSPs in 47 Tuc (diamonds) and field (triangles). 
MSP in NGC 6397 (square) is labeled as is -J in 47 Tuc. Values for MSPs in 
field and M28 are from Table 2, where Lradio = SuW^^- 47 Tuc MSP 5i4oo 
values are from Camilo et al. (2000). Correlation line shown with 47 Tuc-C 
excluded (see text). 

A significant number of additional (>, 15-70) MSPs already 
detected in the Chandra data but not yet detected in the deep 
Parkes radio surveys is quite possible, given the relatively flat 
dependence of soft x-ray luminosity, Lj,, on effective radio lu- 
minosity, Lradio = Suood^, where Sim is the flux density at 
1400 MHz and d the MSP distance. This relation is shown in 
Figure 7 for 47 Tuc, NGC 6397, M28 and field MSPs. The 
47 Tuc MSPs have a distinctly flatter dependence on L,adio> 
wifli logarithmic slope 0.00±0.23 for aU the 47 Tuc MSPs or 
0.25±0.15 with 47 Tuc-C excluded from the fit, as shown in 
Figure 7. The logarithmic slope is 0.99±0.17forthe field MSPs 
(and M28), which is at least partly affected by sample selection 
bias. Thus it appears that in 47 Tuc MSPs that are fainter at 
radio wavelengths are not as much fainter in x-rays and should 
be more readily detected in Chandra data than in Parkes data. 
The ratio of 47 Tuc MSPs Ukely detected in x-rays (~35-90) 



vs. radio (20) can constrain the relative beaming fractions since 
both samples are detected above the same x-ray luminosity limit 
(Lx = 10^*^ erg s"'). If the x-rays are not beamed, the impUed 
radio beaming fraction is ~0.2-0.5. 

We cannot (yet) isolate the candidate MSPs among the unre- 
solved component, but as a possible aid to future MSP identifi- 
cations, we report the positions and fluxes of the 19 Soft sources 
detected which are not identified with possible BY Dra or CV 
candidates. In Table 5 we give the source positions, detected 
counts, and approximate luminosities (again, assuming for 
consistency a kT = 0.2 keV BB spectrum in the 0.5-2.5 keV 
band). These are all possible MSP candidates (though some 
BY Dras and a few CVs are likely included) and thus appro- 
priate for use as trial positions in timing solutions for detected 
radio MSPs. 

4. COMPARISON WITH MSP IN NGC 6397 

We now compare the properties of the 47 Tuc MSPs with 
those derived for the single MSP in NGC 6397. This MSP, 
PSR J1740-5340 (D'Amico et al. 2001a), hereafter 6397-A, 
was recently located by a pulse timing solution (DPMOl) and 
identified with an optical counterpart (Ferraro et al. 2001) pre- 
viously identified with a main sequence binary or BY Dra sys- 
tem (TGEOl). Chandra observations of NGC 6397 (GHEOlb) 
detected 20 sources within 1 ' of the cluster center, including 
the MSP. Because the Chandra sources included 8 CVs firmly 
identified with HST optical counterparts displaying Ha emis- 
sion (GHEOlb), the Chandra-HST boresight could be measured 
to ~0f'01 in RA and Dec. Thus the astrometric identification 
of the Chandra source U12 (GHEOlb) with the apparent BY 
Dra binary WF4-1 found by TGEOl is secure (the Chandra- 
HST positional offsets are given in Table 1, and the Chandra- 
radio offsets are given in GHEOlb). Since this binary has an 
optical photometric period 1.34d (TGEOl) which is consistent 
wifli fliat of 6397-A (D'Amico et al. 2001a), and is included in 
the ~1 " astrometric uncertainty of the HST guide star vs. ra- 
dio timing reference frames, it was identified as the probable 
optical counterpart by Ferraro et al. (2001). In addition to the 
precise Chandra positional identification of the 1.34 d binary, 
we confirm the MSP identification by its x-ray temporal vari- 
ability which was mentioned by GHEOlb as being consistent 
with modulation at the binary period. 

In Figure 8 we show the x-ray light curve of 6397-A (source 
U12) plotted against the binary phase as determined from the 
ephemeris of DPMOl. The count rate increases smoothly by a 
factor of ^2 at phase ^0.4, or just before the MSP comes out 
of radio eclipse (which is centered at phase 0.25; phase is the 
time of ascending node), removing any doubt that the MSP is 
the x-ray counterpart. The fact that the flux is not zero during 
eclipse (though Chandra coverage only begins at phase ^0.1) 
suggests that the x-ray source is extended and larger than the 
companion size, which has radius ^1.3-1.8^0 and nearly fiUs 
its Roche lobe (Burderi et al. 2002). Thus the x-ray emission 
is likely to be from the relativistic wind from the MSP interact- 
ing with the mass loss from the companion. Any thermal (BB) 
component from the MSP itself must be ^0.5 the total flux, 
as constrained by the flux increase out of eclipse. The hard- 
ness ratios during (HR1=4.7![^, HR2=0.42±0.16) and outside 
eclipse (HRl=3^Q g, HR2=0.57±0.2) remain the same within 
the errors, implying the modulation is not due to increased Nh- 
More sensitivity and phase coverage is needed to confirm and 
constrain the modulation. 



9 



The x-ray counts detected in each of the Soft, Medium and 
Hard bands are given in Table 1, together with the E value from 
DPMOI, which assumes that, given its relatively large projected 
offset of 0'6 (^11 Tc) from the cluster center, the cluster accel- 
eration may be ignored. MSP 6397-A is clearly very hard, with 
hardness ratios HRl =3.6 and HR2 = 0.5. We do not plot it in 
Figure 2 since the very different absorption column (A^^^ = 1.0 
X 10^^ cm~^, or >,4x that for 47 Tuc) precludes direct compari- 
son with the models plotted for 47 Tuc, but it is consistent with 
a PL with photon index ^ 1.5 or TB with kT >,5 keV. 

However we can compare the derived E and characteristic 
age, and their relation to L,, for this MSP with the 47 Tuc sam- 
ple. In Figures 3 and 4 we plot 6397-A but do not include it with 
the linear fits which were derived for the 47 Tuc MSPs alone. 
It is clear that 6397-A is consistent with the 47 Tuc sample and 
correlations shown in both figures. Including it on the Lj vs. 
E correlation would give a combined (47 Tuc and NGC 6397) 
Ljc(xE relation with f3 = 0.48±0. 15, even including the errors 
on all the points (despite the bias towards steepening the slope 
in doing so, as mentioned before). This is at first surprising, 
given the very different spectral shape and thus emission region, 
as well as (much) smaller spindown age, log(P/2P) = 8.5, for 
6397-A. However the eclipsing MSP in 47 Tuc with the hardest 
spectrum, 47 Tuc-J, already hints that the x-ray spectral shape is 
not critical to the vs. E correlation, since -J does not depart 
from the correlation line. Thus there seems to be a relatively 
fixed (0.5-2.5keV) for both the 47 Tuc MSPs and 6397-A 
which emerges for a given spindown power E regardless of the 
emission mechanism. It appears that the conversion of E into 
thermal (BB) emission from the NS and emission probably tied 
to the relativistic pulsar wind are both fixed by the total polar 
cap current which, from the models of HM02 is in turn propor- 
tional to E (see §5). 
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Fig. 8. — X-ray count rate (0.5-4.5 keV) vs. orbital phase (from ephemeris 
of DPMOI) for the MSP in NGC 6397, source U12 of GHEOlb, observed for 
49ksec beginning at UT 3 1 .647 July, 2000. The source is persistently eclipsed 
in the radio at phases 0.05-0.45. The x-ray detection at this phase, and increase 
at phase >,0.4, imply at least part of the x-ray emission is extended. 

Finally, we note that the large radial offset of 6397-A in the 
cluster, at radius ^11 r^, is much larger than any of the de- 
tected, or plausible new candidate, MSPs in 47 Tuc. DPMOI 
and GHEOlb noted it was likely ejected from the core, perhaps 
in a second exchange encounter since the MSP companion is 



not the expected He WD or low mass (ablated?) remnant of 
a secondary in an LMXB (as for all previously identified MSP 
companions) but rather is apparently just off the main sequence. 
GHEOlb noted that 6397-A is approximately aligned with the 
central "line" of blue stragglers as well as several other of the 
Chandra sources — suggesting possible asymmetric ejection of 
binaries from the cluster core (possibly induced by cluster rota- 
tion; R. Spurzem, private communication). Since NGC 6397 is 
a post core collapse cluster with a power law cusp measured out 
to ^60" (Lugger, Cohn and Grindlay 1995) that is expected to 
contain scattered binaries (and NSs), the MSP has surely been 
scattered. 

Burderi et al. (2002) show that 6397-A could have evolved 
from an initial ~14h binary period to its present 34 h period, 
near the critical value where pressure from the MSP wind ex- 
ceeds that of matter overflowing the Roche lobe (though inter- 
mittent accretion could resume; cf. §5). However the scatter- 
ing encounter probably occurred within the past ^10^ yr since 
the central relaxation time is only ~10^ yr (Harris 1996) and 
still <,10^' yr at the current offset position. An upper limit to 
the scattering time is ^10^~^y, the likely time since core col- 
lapse and also the likely MSP "age" (Figure 4a). The scattering 
timescale is thus shorter than the binary or nuclear evolution 
timescales leading to the present system, and the interaction 
must have altered the binary parameters. Most probably, the bi- 
nary hardened in the encounter if the present secondary was the 
original donor. This makes it less likely that the secondary, now 
just inside its Roche lobe, could have filled its Roche lobe prior 
to the encounter and spun the MSP up in an LMXB phase. If 
mass transfer occurred only after the hardening induced by the 
encounter, there would not be enough time to spin up the MSP 
given the m ~10""' M© yr"' for the hkely intermittent LMXB 
phases (Burderi et al. 2002) and flie MSP "age". 

If, however, the current secondary was exchanged into the bi- 
nary displacing what was the original donor (probably a much 
lower mass HeWD, as in most other MSPs), there would not 
be the MSP formation-timescale problem. Even more interest- 
ing, this may help explain a puzzle pointed out by Orosz & van 
Kerkwijk (2002) after this paper was submitted: the measured 
(DPMOI) E for 6397-A would suggest significant photometric 
heating of the binary companion (by the MSP wind) though 
no visible heating modulation (only eUipsoidal variation) is ob- 
served. The "excess" E might be a factor of 2-3, depending on 
uncertain heat transport in the secondary. Orosz & van Kerk- 
wijk suggest the large E may be due to a P value dominated by a 
triple companion (easily tested by longer term MSP timing). If, 
instead, the present secondary has been exchanged for the origi- 
nal, and it is marginally Roche lobe filling (Burderi et al. 2002), 
then accretion (and LMXB) episodes could lead to spindown 
episodes of the MSP - at least for those encounters that are ret- 
rograde but possibly also for the very misaligned prograde case. 
Depending on uncertain crust-core coupUng timescales, an in- 
creased P (negative; or positive for the nearly aligned prograde 
case) may persist for the MSP between intermittent accretion 
episodes. This enhanced P in an "old" MSP with low B field 
would not produce the same enhanced pulsar wind (and thus 
heating) as a nominal "young" (higher B field) MSP since P is 
now dominated by internal, rather than electromagnetic, rota- 
tion loss. The true age of 6397-A would thus be correspond- 
ingly older. Of course reduction of E would move 6397-A off 
the apparent correlation lines (e.g. Figs. 3, 4a) defined by the 
47 Tuc MSPs, though if (only) a factor of -^2-3 reduction in 
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E this might be compensated if only the thermal component of 
Ljc (not yet isolated) were plotted. 

A second example of a possible "double-exchange" MSP, 
with a probable main sequence secondary, has been discovered 
with the optical identification of the MSP 47 Tuc-W: Edmonds 
et al. (2002) report (after this paper was submitted) the identi- 
fication of Chandra source W29 (cf. GHEOla) with a V=22.3 
star in the central core ( <,0.2re from the center) of 47 Tuc with 
optical modulation identical in both period (to within 0.5 s) and 
phase (to within 1.2 min) with the 3.2 h binary modulation of 
47 Tuc-W. The x-ray luminosity and hardness ratios for the 
Chandra source (W29) are very similar to those for 6397-A. Al- 
though this is a certain MSP identification, its one-time (^4h) 
radio detection (Camilo et al. 2000) precludes measurement of 
its P and thus E value (possibly enhanced?) so that it cannot 
be included (now) in our analysis of the full radio MSP sam- 
ple in 47 Tuc. Its identification does show, however, that our 
estimates for the total MSP sample in 47 Tuc as derived from 
the soft source distribution are likely to be lower hmits. In- 
deed, Edmonds et al. (2002) also find in the core of 47 Tuc 
a second possible MSP with main sequence companion as the 
optical counterpart of Chandra source W34, with orbital period 
97.45m (not matched to any of the radio MSP sample). 

5. DISCUSSION 

One of the key results of this study is that MSPs in the globu- 
lar clusters 47 Tuc and NGC 6397 may have a less efficient con- 
version of rotational spindown energy (E) into soft x-rays (L;^) 
than most field MSPs, even those with correspondingly low val- 
ues of their magnetic field at the light cylinder. For Lj oc E ^, 
the 47 Tuc-NGC 6397 samples are fit by = 0.5 ±0. 15 whereas 
the field and M28 sample are consistent with the value /3 = I 
found for non-recycled pulsars. The 4 field MSPs with largest 
E and Z?ic (cf. Figs. 3 and 5) show correlations with L, which 
extrapolate to the low E and £ic sample, whereas the correlation 
shown by the 47 Tuc-NGC 6397 MSPs does not extrapolate to 
higher E and Bic systems. It is possible that the difference is just 
due to the 4 MSPs with >,10^^ G being dominated by non- 
thermal beamed emission (e.g. since B1821-24 has no binary 
companion, its hard spectrum is not from a MSP wind interact- 
ing with a companion, as plausibly is the case for 47 Tuc-J and 
6397-A). However, then the large deviations from the 47 Tuc 
correlation of the field MSPs at lower Bic must be explained as 
systematic errors in distance, which is unlikely for at least two 
of them. Perhaps more convincing, the IE vs. r relations 
for the field vs. 47 Tuc MSPs (Fig. 4b) show significant devia- 
tions for the low systems. What could be different about the 
MSPs in these two globulars vs. field MSPs? 

The polar cap (PC) heating models of HM02 can predict the 
observed Lj^ -E ^ relation, with (3 '~0.5, for MSPs with spin- 
down ages T = P/2P ^10** yr that are no longer able to pro- 
duce pairs by curvature radiation (CR) 7-rays (which convert 
to pairs in the pulsar magnetic field) but instead by 7's from 
inverse Compton scattering (ICS) of thermal photons from the 
NS surface. >From the numerical models for PC heating by re- 
turning positrons, with total heating luminosity L+ for the ICS 
regime and the resulting prediction for 'L+IE vs. t (HM02, Fig. 
8b), we find the approximate empirical relations h^lE oc t^P^" 
which yields L+ oc E and L+ oc t"' 

-l/p2(l-a) 

usmg T oc 

l/E f 2. For MSPs with r < 1-2 x 10^ yr (beyond which the 
expected L+/E decreases as the MSP approaches the pair death- 
line), this gives a -^0.4-0.5, or L+ oc £ or /? = 1 - a = 0.5- 



0.6. The model predicts (for f = 5 ms and PC temperature T = 3 
X 10^ K) log(L;t /E ) ~-4.5 -I- 0.41og(T/10'' yr) and so is consis- 
tent with the approximate trend from 6397-A to the apparently 
youngest 47 Tuc MSPs. However, several key differences are 
evident: (i) the 47 Tuc MSPs (and 6397-A) show an unbro- 
ken power law dependence on age, with L^; oc L+ oc r""^ (Fig. 
4a), well beyond the expected steepening (HM02, Fig. 9b) at 
T ~10^ yr, (ii) the IE relation (Fig. 4b) Ukewise continues 
to increase beyond the expected flattening at similar r ~ lO'yr 
(HM02, Fig. 8b). In contrast, the field MSPs in Figures 4a,b 
are more consistent with the HM02 models, having negative or 
flat slopes, respectively, at r >,10^yr. However, the factor of 
^3-lOxexcess of Lj IE over those for the 47 Tuc MSPs and 
the ICS models also suggests these systems have an additional 
source of \-,x : probably non-thermal emission. It is interesting 
that accounting in Figure 4b for the approximate dependence 
on pulsar period, IE oc P^" (derived from a fit to Fig. 8b of 
HM02), makes the field MSPs have a (slightly) more negative 
slope (-0. 16 ± 0.23 vs. -0.06 ± 0.21), as expected from HM02 
for their age, whereas the 47 Tuc MSPs have a sUghtly more 
positive logarithmic slope (0.59±0.10 vs. 0.53±0.13). Finally, 
(iii) the small effective BB emission areas {r^ <,0.14-0.56km 
for dipole polar caps; smaller still (each) for multipole caps) 
are significantly smaller than the effective radius of the annular 
polar cap emission region predicted (HM02 Fig. 7). 

We suggest that these differences may point to altered sur- 
face magnetic field configurations for the 47 Tuc (and 6397-A) 
MSPs to enable such efficient PC heating at large age. If their 
surface magnetic fields have significant multipole components, 
the effective B-field curvature is greater and the pair production 
front moves closer to the NS yielding more efficient positron 
return to heat the polar caps. The magnetic field footprint on 
the NS becomes smaller, giving the small emission areas. If 
the non-thermal radio emission is produced, as in most models, 
from near the light cylinder where the field is dominated by the 
dipole component, then its properties are not significantly al- 
tered and the radio pulse profiles are otherwise similar to field 
MSPs (as observed; Camilo et al. 2000). Clearly this multipole 
field hypothesis must be tested with models (only dipole fields 
have been considered for pulsar emission and PC heating — 
e.g. HM02). 

Another possibility is that the NSs in the 47 Tuc MSPs (and 
6397-A) are systematically more massive or compact than those 
in the field MSPs. The models of Harding, Muslimov, & Zhang 
(2002) show that the ICS regime becomes longer lived (extends 
to lower P values) for more massive or more compact NSs. Al- 
though our radial distribution results (§3.2) imply MSP masses 
(total) in the range ~ 1.1-1. 4 Mq , suggesting the more com- 
pact NS interpretation may be preferred, more massive NSs are 
allowed if in fact the MSP distribution is "puffed up" by binary 
scattering. 

The hard spectra in a few of these "thermal" MSPs (47 Tuc- 
J and 6397-A) must have a different origin. In 6397-A, the 
possible modulation phased with the radio eclipse of the hard 
source suggests the hard source is due to shocked gas hfted 
from the binary companion by the MSP wind pressure (Bur- 
den et al. 2002). If the hard spectrum is a power law (thermal 
bremsstrahlung is ruled out by the (5DM constraint [§2.1] un- 
less the emission is strongly aspherical), and the shocked wind 
density and thus emission measure is proportional to the under- 
lying MSP thermal luminosity, then the approximate agreement 
of these hard spectrum, but low \-,x, systems with the (3 = 0.5 
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law exhibited by the soft (thermal) MSPs would imply scaling 
of the shocked gas emission with the thermal (NS) emission. 
As mentioned above, this might be expected if both the MSP 
wind and PC heating are proportional to the PC current. 

What is special about globular cluster MSPs to drive them 
towards altered fisurf (and thus fit) values or configurations 
or possibly larger NS mass (or compactness)? Unlike MSPs 
in the field, those in dense cluster cores have a possibiUty 
of being driven back into contact and an accretion phase, as 
a re-cycled LMXB (from a MSP)! Renewed accretion (and 
the 6397-A companion is presently close to filhng its Roche 
lobe; cf. DPMOl) would Ukely continue the B-field burial pro- 
cess thought to be responsible (e.g. Romani 1990) for field de- 
cay from the > 10" G fields at NS birth to the < val- 
ues typical of MSPs. This could lead, in turn, to altered B- 
field configuration (multipole field components), particularly 
for a secondary re-exchanged into an MSP (with random spin- 
orbital encounter angular momentum). It would also systemat- 
ically increase the NS mass beyond that for field MSPs. The 
47 Tuc MSPs have spent their lifetime r > 1 Gyr in a dense 
(n ~ lO^pc"^) cluster core, where they undergo scattering in- 
teractions with both single stars and other binaries. Such scat- 
tering causes the binding energy x of a hard binary to secu- 
larly increase ("hard binaries harden") at a rate x ~ \ma^ np^a 
(Goodman & Hut 1993), where m is the single star mass, a is 
the one-dimensional velocity dispersion, and po = Gm/cP' is the 
impact parameter for a 90° deflection (~ 10 AU in the core of 
47 Tuc). Thus, the binding energy of a typical hard binary in- 
creases by of order ma^ per core relaxation time (~ 10** yr for 
47 Tuc). Over the > 1 Gyr Ufetime of the MSP binary system, 
the total relative binding energy change is xt/x ^ 1%. Thus, 
while angular momentum transfer to the secondary acts to de- 
tach the binary, scattering events tend to drive it back toward 
contact. 

A complication in this picture is that since MSPs are ex- 
tremely hard binaries (x ^ 1 000 mcr^), x is adiabatically in- 
variant in weak scattering events where the distance of closest 
approach much exceeds the binary separation (e.g. Heggie & 
Hut 1993). Thus, the secular increase in x is largely the re- 
sult of infrequent strong scattering events, which are likely to 
result in an ejection of the binary from the cluster core (e.g. 
6397-A) or even the cluster. Other unresolved issues include 
the effects of binary-binary scattering, which probably "bump" 
the MSPs back into contact more effectively than does binary- 
single scattering, and orbital eccentricity evolution during scat- 
tering events, which will also tend to drive some binaries to- 
wards contact. Nevertheless, it appears plausible that a typical 
old (r > 1 Gyr) MSP in a dense cluster core might undergo 
one or more MSP-LMXB recurrence cycles. Thus a few per- 
cent of the MSPs at any one time may be in (or near) this re- 
current LMXB phase. This might explain the puzzling lumi- 
nous qLMXBs X5 and X7 (GHEOla and Heinke et al. 2002) in 
47 Tuc: they could be recently (~ lO'' "^yr) revived qLMXBs, 
in which their underlying MSP nature is hidden by a heated 
(by the MSP wind) outflow from the disk and Roche lobe of 
the probable main sequence secondaries (at least for X5, which 
may then be similar to 6397-A and 47 Tuc-W; cf. §4). 

The (possibly) younger system in NGC 6397 has the advan- 
tages of having been recently scattered out of a still higher 
density (~10^ pc~^) core collapsed cluster core and possibly 
having exchanged its companion. Either or both would Ukely 
have restored an accretion phase. Thus the MSP need not be 



just "born," as suggested by Ferraro et al. (2001); it may in- 
stead have just been reborn. In contrast, the MSP in M28 is 
both ~10x younger, single and in a lower density core and so 
is unlikely to have gone through a renewed accretion phase. 
This MSP-LMXB recycling scenario would imply: (i) a pop- 
ulation of MSPs with main sequence or red straggler secon- 
daries (either primordial or re-exchanged) like 6397-A which 
may be largely ecUpsed (in radio) by mass loss and dominated 
by hard shocked gas emission, (ii) a strong dependence for this 
multiply-recycled population to occur primarily in clusters with 
largest central mass density po = Ho m and core radius r^, and 
thus total interaction rate (F oc pl'^r^) such as Terzan 5, and 
(iii) that multiple MSP-LMXB phases may be more likely, thus 
yielding Bsud burial and altered Bic and suppressed magneto- 
spheric emission, in MSPs with main sequence companions. 
MSPs with either main sequence or slightly evolved secon- 
daries (e.g. 47 Tuc-W or 6397-A; cf . § 4) are also most likely to 
contain re-exchanged secondaries which may have anomalous 
P values (as might 6397-A) from the braking torque supplied 
by a mis-aligned crust-core spin coupling due to the renewed 
accretion phase(s) from the re-exchanged secondary. Even an 
original secondary, if driven back into contact or if marginally 
filling its Roche lobe, could give rise to enhanced P values: 
when accretion stops, the crust is spinning faster than the core, 
which then exerts a braking torque on the crust (and B-field) 
of the MSP. The timescale for equilibrium to be established is 
uncertain but likely to be relevant for MSPs with companions 
closest to filling their Roche lobes (hence the main sequence 
or subgiant systems) which have most recently accreted. These 
main sequence systems, which likely evolve into the short pe- 
riod (~1.5-5.5h) group identified by Camilo et al. (2000), also 
have longer lifetimes to undergo binary scattering than those 
spun-up from sub-giant companions and now in the longer pe- 
riod (~0.4-2.4d) group with He- WD companions. However 
this Ufetime is partly offset by the larger binary cross sections 
for the He-WD group, so both are expected to undergo scatter- 
ing in dense cores. The generally larger orbital eccentricities 
for He-WD secondary MSPs in globulars vs. those in the disk 
has previously been noted (Phinney 1992) as evidence for inter- 
action of these MSPs with cluster stars. The large eccentricity 
of 47 Tuc-H (FCLOl) suggests it has undergone such an en- 
counter. 

6. CONCLUSIONS 

In this paper we have examined the x-ray properties of the 
full MSP populations in 47 Tuc and NGC 6397 localized via 
radio timing. The large sample of MSPs resolved by Chandra 
in 47 Tuc, and the single system in NGC 6397, allow a number 
of conclusions to be drawn: 

• The MSPs in 47 Tuc and NGC 6397 are surprisingly 
uniform in despite a much larger variation in E. The 
hx-E correlation appears to be significantly flatter than 
for a (smaller) sample of MSPs in the field as well as the 
isolated (and young) MSP in M28. The correspondingly 
flatter l^x-Bic correlation suggests that the 47 Tuc (and 
NGC 6397) MSPs may have B,^^ and thus configu- 
rations systematically different (e.g. multipole compo- 
nents) from field MSPs and which prolong thermal and 
suppress magnetospheric x-ray emission. 

• The MSPs in 47 Tuc are dominated by thermal emis- 
sion sources, consistent with BB emission from the pul- 



sar polar caps. At least one (47 Tuc-J) must also have 
significant non-thermal emission, which may arise from 
the pulsar wind interacting with the companion (see be- 
low). 

The radial distribution of MSPs in 47 Tuc is consistent 
with the predictions of a multi-mass King model for a 
distribution of MSPs with mass ratio ^ 2 times the 
typical cluster stellar mass in the core, i.e. a MSP sys- 
tem mass of 1 .4 Mq . 

From the distribution of Chandra sources with similarly 

soft x-ray colors, as well as the total counts detected 
in Soft and Medium bands vs. cluster radius, the num- 
ber of MSP candidates with ^1 x 10""' erg s^' is es- 
timated to be ~35-90, allowing for the plausible con- 
tamination of BY Dra and soft CV systems and incom- 
pleteness of the detection for faint sources near the clus- 
ter center. The relative number of MSPs detected in 
the complete soft x-ray vs. radio samples suggests the 
radio emission is beamed with beaming factor ~0.2- 
0.5 (times any soft x-ray beaming, which is presumably 
isotropic). 

A larger population of still fainter MSPs (Lx < 10^0 
erg s~') could be detected in x-rays than in the radio, 
given the relatively flat Li-Lradio correlation seen for 
47 Tuc MSPs. If the excess of Soft vs. Medium band 
counts in the cluster core is due to fainter MSPs, they 
must have relatively lower luminosity and lower tem- 
peratures than the identified MSPs. 

The one (binary) MSP system identified in the core- 
collapsed globular NGC 6397 is consistent with the L^- 
E relation of the 47 Tuc sample despite its very differ- 
ent spectrum, which must be dominated by non-thermal 
or bremsstrahlung emission. Its harder x-ray emission 
must (largely) arise from an extended source, possibly 
resulting from interaction of the relativistic wind from 
the MSP with mass loss from its companion star. The 
fact that it is still consistent with the Li oc £ '/^ rela- 
tion found for the thermal MSPs in 47 Tuc suggests that 
both the pulsar wind luminosity and polar cap heating 
are governed by the total polar cap current (oc £ cf. 
HM02). If P is enhanced by renewed accretion or sec- 
ond exchanges, the electromagnetic spindown rate and 
corresponding reduced E may correlate with the thermal 
component of (to be measured) as for 47 Tuc. 

The enhanced thermal luminosity, small emission areas 
and extension to larger MSP ages, for the 47 Tuc MSPs 
vs. the HM02 models suggests the surface magnetic 
field might be multipole (dipole at the light cylinder ra- 
dius and probable radio emission region), possibly by 
field perturbations from repeated accretion events due to 
encounters or re-exchanges of the secondaries of these 
MSPs in dense cluster cores. Alternatively, the extended 
pair heating (and thermal x-rays) could arise (Harding, 
Mushmov, & Zhang 2002) if the NSs in these MSPs are 
more massive or (given their radial distribution) more 
compact. If re-exchanged secondaries are present (as 
perhaps for 6397-A and 47 Tuc-W), apparent P and thus 
E values may be enhanced by core-crust coupling drag 



effects over their purely electromagnetic spindown val- 
ues which give rise to heating their companions. Some 
evidence for a similar "excess" E value is available for 
47 Tuc-W (Edmonds et al. 2002), but as with 6397-A 
these estimates are subject to uncertain heating efficien- 
cies and transport effects (which might heat the backside 
of the tidally-locked secondary). 

• NGC 6397 appears to be deficient in its total MSP pop- 
ulation per unit mass compared with 47 Tuc (GHEOlb), 
if most cluster MSPs are like the 47 Tuc sample and 
dominated by soft thermal sources. Only 3-5 such ob- 
jects are found in the sample of 20 Chandra sources 
detected, with L^: limits a factor of 2 lower than in 
47 Tuc (GHEOlb). Most of the other Chandra sources 
in NGC 6397 with comparable x-ray colors (or actual 
spectra) are identified with CVs (optically; cf. GHEOlb 
and references therein) and thus not likely to be addi- 
tional MSPs. However it is possible that a few of the 
main sequence binaries, or BY Dra candidates (TGEOl) 
could also harbor MSP primaries as is the case for 6397- 
A. 

• Conversely, the very large possible BY Dra population 
in 47 Tuc (up to ~10^ are possible; R. GilUland, pri- 
vate communication) may contain a significant sample 
of MSPs as may have occurred in 6397-A: those that 
have likely had a cluster main sequence star exchanged 
into the binary, displacing the probable degenerate sec- 
ondary which had evolved from the original mass donor 
to spin up the MSP. If the 6397-A system with such a 
secondary (now evolved or bloated; cf. § 4) is typi- 
cal, and mass loss is more readily driven off its near 
Roche-lobe filling secondary by the MSP wind pres- 
sure (Burderi et al. 2002), then these systems will 
likely be dominated by shocked gas and appear hard. 
Since these main sequence secondary MSPs will typi- 
cally have larger MSP heating (given their larger secon- 
daries) and thus mass loss, they may be largely eclipsed 
until their secondaries have been ablated to become the 
very low mass group (^O.O2M0) found (Camilo et al. 
2000) in the 47 Tuc MSPs. Thus the total MSP sample 
in 47 Tuc would include faint hard sources as weU as 
the soft, thermal-dominated, systems and the total MSP 
population would exceed the ~35-90 estimated from the 
soft source population alone. 

We shall obtain a very deep (300ksec) cycle 3 Chandra ob- 
servation of 47 Tuc to probe the complete MSP population as 
well as to better constrain and measure x-ray colors and actual 
spectra to test emission models. The hard spectra MSPs (e.g. 
47 Tuc-J) may be variable (as is 6397-A), and hard components 
should be detectable in other eclipsing MSPs (e.g. 47 Tuc-O). 
The interlocking puzzles of compact binary populations (e.g. 
MSPs vs. BY Dra systems) and compact object populations 
(MSPs vs. CVs, as measures of NS vs. WD populations) in 
globular clusters offer new insights into the stellar and dynam- 
ical evolution of these oldest stellar systems in the Galaxy. 
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Table 1 



Chandra Results on MSPs in 47 Tuc and NGC 6397 



MSP 


ARA^ 


ADec^ 


Softcts'' 


Mediumcts'' 


Hardcts'' 


log(LJ'= 


log(£)'' 




(arcsec) 


(arcsec) 


(0.2-1 keV) 


(1-2 keV) 


(2-8 keV) 


(0.5-2.5 keV) 




Al Tun C 






1 
1 


n 
u 


n 
u 






47 Tuc-D 


-0.17 


0.19 


5 


7 


1 


30.3 


33.7 


47 Tuc-E 


0.01 


-0.07 


11 


13 





30.6 


34.5 


47 Tuc-F 


0.32 


0.06 


10 


9 


1 


30.5 


34.4 


47 Tuc-G 


0.05 


0.03 


4 


4 





30.1 


33.7 


47 Tuc-H 


0.03 


0.01 


5 


4 





30.1 


33.2 


47 Tuc-I 


0.16 


0.01 


5 


5 


1 


30.2 


33.9 


47 Tuc-J 


-0.03 


-0.07 


3 


5 


5 


30.3 


34.4 


47 Tuc-L 






12(10) 


5 


6(0) 


30.4 


34.2 


47 Tuc-M 


0.36 


-0.01 


2 


3 





30.1 


33.5 


47 Tuc-N 


0.01 


0.10 


5 


6 





30.2 


34.1 


47 Tuc-0 


0.10 


0.19 


17 


11 





30.6 


34.4 


47 Tuc-Q 






5 


3 





30.1 


34.3 


47 Tuc-S 


-0.26 


-0.39 


5 


4 


1 


30.2 




47 Tuc-T 






4 


3 





30.0 


33.8 


47 Tuc-U 


0.04 


-0.06 


4 


7 


1 


30.3 


34.2 


6397-A 


-0.05 


0.09 


11 


40 


20 


30.9 


35.1 



^ ARA and ADec offsets are Chandra wavdetect minus precise MSP radio (timing) position (47 Tuc) 
and Chandra WAVDETECT minus HST position (NGC 6397). The close pair 47 Tuc-F, -S are marginally 
resolved (see text) and thus have larger offsets; offsets are not listed for -L, -Q, -T and -C, without 
WAVDETECT positions duc to crowding or limited counts (see Fig. 1). 

''Counts in each band derived from 1 " aperture on radio MSP position. MSPs 47 Tuc-G and -I are 
indistinguishable; counts are tentatively divided based on respective E values. MSPs 47 Tuc-F and -S 
detected combined with wavdetect; counts estimated separately (see text). 47 Tuc-L suffers severe 
crowding; values in parentheses indicate best estimates removing other source contamination. 

•^L^ values (in erg s"') are computed (PIMMS) for an assumed blackbody spectrum with kT=0.22keV 
(see text) and include exposure correction factors of 2.5, 1 .4 and 1 .9 for 47 Tuc-C, -J, and -M, respectively 
(all other MSPs detected with uniform exposure factor 1.0). Approximate errors on log Ljc are ±0.3 for 
47 Tuc-C and ±0.2 for the others. 

values (in erg s~') are derived for the 47 Tuc MSPs assuming the MSPs are at 3D radii appropriate 
to their dT>M, determined using a constant gas density in the cluster (Freire et al. 2001b), and a King 
model for the cluster potential with central velocity dispersion 11.6kms"'. E value for 47 Tuc-S is 
indeterminate. E for the MSP 6397-A, at relatively large radial offset in NGC 6397, is derived assuming 
the cluster acceleration contribution is small. 
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Table 2 

Parameters for MSPs Previously Detected in X-rays 



MSP 


P 


Pi 


d 


>5l400 


F;c^ References 




(ms) 


(xlO-2") 


(kpc) 


(mJy) 


(xlO""ergs-i cm-2) 



J0030+0451 


3.05 


<1.0 


0.23 


0.6 


2.5 


1,2 




J02 18+4232 


2.32 


7.78 


>5.9 


0.9 


1.5 


3,4 




J0437-4715 


5.76 


1.86 


0.14'' 


137 


19. 


4,5, 


6 


J0751+1807 


3.48 


0.80 


2.02 


1.0 


0.83 


4,7 




J1012+5307 


5.26 


0.97 


0.84 


2.8 


0.49 


4, 7, 


8 


J 1024-07 19 


5.16 


<0.30 


<0.20'' 


0.9 


0.2 


4, 6, 


9 


J1744-1134 


4.07 


0.71 




2.0 


0.2 


4,6, 


10 


B1937+21 


1.56 


10.6 


>3.6'' 


16 


3.7 


7,9, 


11 


B 1957+20 


1.61 


1.15 


1.53 


0.35 


3.0 


4, 7, 


9 


J2124-3358 


4.93 


1.30 


0.25 


2.6 


2.5 


4, 6, 


9 


B1821-24 (M28) 


3.05 


162. 


5.5" 


1.9 


5.2 


4,7 





References. — 1 (Lommen et al. 2000); 2 (Becker et al. 2000); 3 (Navarro et al. 1995); 4 
(BT99); 5 (van Straten et al. 2001); 6 (Toscano et al. 1998); 7 (Taylor et al. 1995); 8 (Lange et 
al. 2001); 9 (Toscano et al. 1999b); 10 (Toscano et al. 1999a); 11 (Takahashi et al. 2001). 

''X-ray flux in ROSAT (0.1-2.4keV) band. Value forB1937+21 in ASCA (2-lOkeV) band. 

''Accurate distances/limits. Other estimates are obtained primarily from dispersion measures 
together with the Taylor & Cordes (1993) electron density model, and are rather uncertain. 



Table 3 

Results from KS Two-Sample Comparisons 





MSP 


Medium 
(raw) 


Medium 
(corr) 


Soft 


0.66 


0.0034 


0.11 


Medium (raw) 


0.16 






Medium (corr) 


0.82 







Note. — The entries give the probabiUty of the null hy- 
pothesis that the two samples are drawn from the same 
underlying distribution. Values of less than 0.05 indicate 
significant differences between the samples. 
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Table 4 

Results from Generalized King Model Fits 



Sample " 


A'' 


rc Cf 


a 


Soft 


44 


15.2± 


3.5 


-3.6 ±0.8 


Medium (raw) 


41 


5.6± 


3.9 


-1.9±0.3 


Medium (corr) 


33 


9.0± 


6.9 


-2.5 ±1.0 


MSP 


16 


13. 1± 


6.9 


-2.4±1.1 


CV 


19 


10.3 ±10.0 


-3.1 ±3.7 


BY + RS 


14 


15.8± 


4.5 


-5.6±4.1 


No ID 


22 


17.5 ± 


3.6 


-3.1 ±0.7 



"Samples were delimited by a radial cutoff of 
250 " except for the MSP group for which 1 00 " was 
used. 

''Core radius for x-ray source population (vs. rc 
= 24.0" for the optical; see text). 

'^BY + RS sample is the combination of BY Dra 
(main sequence) binaries and some possible RS 
CVn (sub-giant) systems. 



Table 5 



Candidate MSPs from Soft Sources 



RA'' 


Dec" 


Medcts 


log(L,)'' 


(J2000) 


(J2000) 


(0.5-4.5 keV) 


(0.5-2.5 keV) 


00:24:13.45(1) 


-72:04:51.27(8) 


33.5 


30.7 


00:24:12.53(1) 


-72:04:41.13(5) 


8.3 


30.1 


00:24:11.64(3) 


-72:05:03.72(8) 


9.1 


30.1 


00:24:07.38(2) 


-72:04:49.41(5) 


41.8 


30.8 


00:24:06.05(2) 


-72:05:01.7(1) 


12.7 


30.3 


00:24:04.64(2) 


-72:04:34.5(1) 


18.7 


30.5 


00:24:03.22(2) 


-72:04:29.74(8) 


15.7 


30.4 


00:24:03.14(2) 


-72:04:26.55(9) 


16.7 


30.4 


00:24:02.63(3) 


-72:05:38.5(1) 


13.8 


30.3 


00:24:02.52(2) 


-72:04:44.1(1) 


20.4 


30.5 


00:23:59.91(2) 


-72:05:04.7(1) 


21.2 


30.5 


00:23:59.45(2) 


-72:04:47.9(2) 


9.4 


30.2 


00:24:01.45(3) 


-72:04:41.6(1) 


9.7 


30.2 


00:24:01.04(2) 


-72:04:13.9(1) 


12.5 


30.3 


00:23:59.45(3) 


-72:04:38.5(1) 


9.9 


30.2 


00:24:05.03(4) 


-72:05:02.8(1) 


9.0 


30.1 


00:24:17.20(2) 


-72:03:03.2(1) 


8.7 


30.1 


00:24:11.50(1) 


-72:02:24.7(1) 


20.6 


30.5 


00:24:42.49(4) 


-72:06:20.6(2) 


15.1 


30.4 



Note. — The first 5 source entries, and the 4th from the bottom of 
the table, are all particularly promising MSP candidates as they are all 
included in the Albrow et al. (2001) HST field and yet lack evidence for 
optical counterparts with either CV or BY Dra characteristics. 

"The RA and Dec have been corrected to the MSP coordinate frame. 



^'Lx values calculated as for Table 1. 



